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Abstract 

The existence of stationary solutions to the Einstein- Vlasov system 
which are axially symmetric and have non-zero total angular momen- 
tum is shown. This provides mathematical models for rotating, general 
relativistic and asymptotically flat non-vacuum spacetimes. If angu- 
lar momentum is allowed to be non-zero, the system of equations to 
solve contains one semilincar elliptic equation which is singular on the 
axis of rotation. This can be handled very efficiently by recasting the 
equation as one for an axisymmetric unknown on R 5 . 
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1 Introduction 



The geometric features of general relativistic and asymptotically flat space- 
times strongly depend on whether the spacetime has non-trivial total angular 
momentum or not. A well known point in case is the Kerr family, a two 
parameter family of stationary vacuum spacetimes which contain a black 
hole. The two parameters are the ADM mass A4 > and the total angular 
momentum C > 0. If C = one obtains the Schwarzschild spacetime, i.e., a 
static, spherically symmetric black hole of mass Ai. Compared to this the 
case C > exhibits a vastly more complicated geometry, and we refer to 
[15| for details. However, most astrophysical objects are not exactly spher- 
ically symmetric, and many rotate about some axis and have non-trivial 
total angular momentum. Hence the mathematical difficulties entailed by 
giving up spherical symmetry and allowing for non-zero angular momentum 
have to be overcome in order to get closer to physically meaningful models. 
There are at the moment only two papers where the existence of rotating 
equilibrium configurations of self-gravitating matter distributions is shown 
in the framework of General Relativity: these are [TO] and [2], where matter 
is modeled as an ideal fluid and as an elastic body, respectively. 

In the present paper we consider matter described as a collisionless gas. 
In astrophysics, this model is used to analyze galaxies or globular clusters 
where the stars play the role of the gas particles and collisions among these 
are sufficiently rare to be neglected. The particles only interact by the 
gravitational field which the ensemble creates collectively, and the general 
relativistic description of such an ensemble is given by the Einstein- Vlasov 
system. The existence of spherically symmetric steady states to this system 
has for example been shown in [14]. In [4] the present authors proved the ex- 
istence of static, axially symmetric solutions which are no longer spherically 
symmetric, but which still have zero angular momentum. In the present 
paper we also remove the latter restriction, a task which, in view of what 
was said above, is not trivial. 

We shall formulate the Einstein- Vlasov system in standard axial coor- 
dinates f 6 K, p £ [0, oof, z £ M, (p £ [0, 2ir]. Following [5], we write the 
metric in the form 

ds 2 = _ c 2 e 2u/c^ dt 2 + e 2» df) 2 + ^2 + (( ^ _ ^2 (L1) 

for functions B,p,u depending on p and z. The reason for keeping the 
speed of light c as a parameter in the metric will become clear shortly. The 
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metric is to be asymptotically flat in the sense that the boundary values 



lim (\u\ + \n\ + \u\ + \B-l\)(j>,z)=0 

\[p,Z)\-HX> 



(1.2) 



are attained at spatial infinity with certain rates which are specified later. 
In addition we require that the metric is locally flat at the axis of symmetry: 



The solutions obtained in [3J are static and have zero total angular mo- 
mentum. In terms of the metric above this means that uj = 0. The quan- 
tity oj is the angular velocity with respect to infinity of the invariantly de- 
fined zero angular momentum observers with worldlines perpendicular to 
the {t = const.} hyper surf aces; see [5]. For a rotating configuration, oj must 
not vanish identically. As in [5] the solutions to the Einstein- Vlasov system 
are obtained by perturbing off a spherically symmetric Newtonian steady 
state using two parameters, 7 = 1/c 2 to turn on general relativity, and A to 
turn on angular momentum. In order to apply the implicit function theorem 
and to make certain solution operators well-defined, it becomes essential to 
handle the linearized w-equation 



for a suitable class of right-hand sides q. The difficulty with this fairly in- 
nocent looking elliptic equation is that the coefficient 3/p blows up on the 
axis of symmetry, where the full solution must remain smooth. It is techni- 
cally very demanding to make this equation as it stands fit into the general 
framework of our approach, cf. the corresponding remark in the appendix. 
However, this equation is nothing but the Poisson equation on M 5 where both 
lo and q are axially symmetric, i.e., they depend on p = \ (xi, X2, x%, x^)\ and 
z = X5. This observation turns out to make the inclusion of non-trivial total 
angular momentum quite neat. We are not aware of a physical background 
for this fact, nor are we aware that this observation has previously been ex- 
ploited in the area of mathematical relativity. It turns out that an analogous 
observation applies to the linearized equation for B which can be turned into 
the Poisson equation on R 4 . This simplifies the proof and improves the result 
also for vanishing angular momentum, when compared with [JJ. One should 
realize that the generalization to non-trivial angular momentum means that 
one moves to a geometrically more complex spacetime. To appreciate the 
fact that the resulting complications are of a genuinely relativistic, geometric 



i/(0, z)/c 2 +n(0,z) = In B(0,z), z£K. 



(1.3) 



dppU + 8 zz uj + - dpUJ = q 



(1.4) 
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nature one should notice that in [13], where an analogous strategy was used 
to obtain axially symmetric steady states in the Newtonian case, i.e., for the 
Vlasov-Poisson system, one and the same proof gives static solutions with 
zero total angular momentum and stationary ones which rotate, depending 
only on which particular ansatz function is chosen. 

Let us now give a formulation of the Einstein- Vlasov system. In a kinetic 
model like the Vlasov equation the particle ensemble is described by its 
distribution function / > which is defined on the tangent bundle TM 
of the spacetime manifold M. Let g a p denote the Lorentz metric on the 
spacetime and let p a denote the canonical momentum coordinates which 
correspond to the chosen coordinates x a on M. The Einstein field equations 

G ap = 8TTc- A T a p (1.5) 

are then coupled to the Vlasov equation 

A»/-r?/^ p0 / = o (i.6) 

via the definition of the energy momentum tensor 

rp -I, ,1/2 f f dp dp 1 dp 2 dp 3 

T a = c \g\ 1 PaPpf • (1.7) 

J m 

Here are the Christoffel symbols induced by the metric, \g\ denotes the 
modulus of its determinant, and m > is the rest mass of the particle with 
phase space coordinates (x a ,pP). The characteristic system of the Vlasov 
equation (jl.6p are the geodesic equations written as a first order system on 
TM. For physical reasons we must require that all particles move forward 
in time, i.e., p a is a timelike, future pointing vector on the support of /. 
Moreover, we make the standard assumption that all particles have the 
same rest mass which we normalize to unity. The distribution function / is 
then supported on the mass shell 

PM = {g a pV a vP = -c 2 m 2 = -c 2 and p a is future pointing} C TM. (1.8) 

It is now important to realize that due to the presence of u in the metric, 
i.e., due to the fact that we want to allow for non-trivial angular momentum 
of the spacetime, the mass shell condition can in general not be used to 
express p° by the remaining variables on TM. It turns out that this can be 
done if and only if the Killing vector d/dt which corresponds to the time 
translation symmetry is timelike everywhere, i.e., 

- g{d/8t, 8/dt) = c 2 e 2u / c2 - p 2 BWe~ 2u/c2 > 0. (1.9) 
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For the solutions which we construct we do a priori not know whether this 
property holds or not. We can if we wish make sure that it does hold so 
that there is no ergosphere. The question whether among the solutions we 
construct there are solutions that do have an ergosphere is open. The vector 
field d/dt + cud/dip is always timelike and can therefore be used to fix the 
time orientation of the spacetime. For the solutions we construct, 

- g(d/dt + ojd/d(p,p a ) = c 2 e 2v / c y > (1.10) 

on the support of / so that all particle trajectories have future pointing 
tangent vectors as desired. We refer to [3] for more background on the 
Einstein- Vlasov system and state our main result. 

Theorem 1.1 There exist stationary solutions of the Einstein- Vlasov sys- 
tem li 1.5\) . Iil.6}) . J 1. 7| ) with c = 1 such that the metric is of the form U.l\) 
and satisfies the boundary conditions Al.fy) , U.3\) . and where the total angu- 
lar momentum is non-zero. 

For the proof of this result the following observation is important. The 
symmetries of the metric imply that the quantities 

E := -g(d/dt,p a ) = c 2 e 2u/c2 p° + p 2 B 2 ue- 2v l c \p z - up ), 
L := g{d/dp,p a ) = p 2 B 2 e~ 2u/c \p 3 -up ), 

are constant along geodesies; notice that 

E = cV^V + uL. 

Here E can be thought of as a local or particle energy and L is the angular 
momentum of a particle with respect to the axis of symmetry. The require- 
ment that p a be future pointing implies that E > on the support of /, 
provided (|1.9p holds, i.e., provided there is no ergosphere. Up to regularity 
issues a distribution function / satisfies the Vlasov equation if and only if 
it is constant along geodesies. Hence we make the ansatz 

f = *(E,L)8(m-l), 

and the Vlasov equation (|1.6|) holds. The 5 distribution on the right hand 
side is to restrict / to particles with rest mass equal to unity; notice that 
the rest mass m is conserved along geodesies as well. If we insert this 
ansatz into the definition (II. 7p of the energy momentum tensor the latter 



5 



becomes a functional T Q( g = T a p(y, B, fj,, uj) of the yet unknown metric func- 
tions v,B,fj,,cj, and we are left with the problem of solving the field Einstein 
equations (|1.5p with this right hand side. 

We will obtain the solutions by perturbing off spherically symmetric 
steady states of the Vlasov-Poisson system via the implicit function theo- 
rem; the latter system arises as the Newtonian limit of the Einstein- Vlasov 
system. We will specify conditions on the ansatz function <3? above such that 
a two parameter family of axially symmetric solutions of the Einstein- Vlasov 
system passes through the corresponding spherically symmetric, Newtonian 
steady state. The parameter 7 = 1/c 2 turns on general relativity and a 
second parameter A turns on the dependence on L and hence axial symme- 
try; L is not invariant under arbitrary rotations about the origin, so if / 
depends on L the solution is not spherically symmetric. Moreover, suitable 
assumptions on the ansatz function $ will force uj to be non-trivial so that 
the solution rotates about the axis p = and has non-zero total angular mo- 
mentum. One should be careful to notice here that there are always particles 
with non-zero angular momentum L provided / is non-trivial and smooth 
on the mass shell, but in general this does not imply that the total angular 
momentum of the whole spacetime is non-trivial. The scaling symmetry of 
the Einstein- Vlasov system can then be used to obtain the desired solutions 
for the physically correct value of c, and not only for large c. 

The idea of employing the implicit function theorem to obtain new equi- 
librium configurations of self-gravitating matter distributions from known 
ones can be traced back to L. Lichtenstein, who investigated the existence 
of non-relativistic, axially symmetric, stationary, self-gravitating fluid balls 
|11[ I12j . His arguments were put into a rigorous and modern framework 
in [9] and extended to the general relativistic set-up in [10] . As mentioned 
above, rotating elastic bodies were considered in [2]. Our approach signifi- 
cantly differs from [21 [10] not only in the matter model, but also in that we 
use the explicit form of the metric stated in ([l.ip . together with a reduced 
version of the Einstein field equations. 

The outline of the paper is as follows. The detailed formulation of our 
main result and the set-up for the application of the implicit function the- 
orem are stated in the next section. In Section [3] we then give a detailed 
outline of its proof. The proof consists of several steps and some of them 
are more or less identical to the corresponding steps in [3] and need not be 
repeated. However, the logical structure of the present proof will be given 
in full detail. In Section [3] we collect some properties of the matter terms 
which will be needed throughout. Section [5] contains information on certain 
Newton potentials which is then used to show that the operator to which we 
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apply the implicit function theorem is well defined. Section [6] explains how 
a solution of the reduced field equations leads to a solution of all the field 
equations. In an appendix we collect a few general results on the regularity 
of axially symmetric functions, and we comment on solving the equation 
(|1.4|) without resorting to the device of moving it into a higher dimension. 

2 Basic set-up and the precise result 

We introduce the parameter 7 = 1/c 2 € [0,oo[. In order to handle the mass 
shell condition effectively it is useful to introduce new momentum variables 

v ° = eVpO, v 1 = e V, v 2 = eV, v 3 = pBe'^ip 3 - ojp°). (2.1) 

This turns the mass shell condition for general m into 

2 2 2/ 0\2 1 1 1\2 1 1 2\2 1 / 3\2 / 0\2 2 , | |2 

—cm = — c [v ) + [v ) + (v ) + [v ) or (v ) = m + "i\v\ 

where v = (v 1 ,v 2 ,v 3 ) G R 3 and \v\ is the Euclidean norm on M 3 . We 
eliminate v° by choosing the positive root which makes sure that (jl.lOp 
holds, i.e., all particles move forward in time. With m = 1 we find that 

E = c 2 e 7 Vl + 7M 2 + "L, L = pBe^v 3 . 

In particular 

E > c 2 e 7 Vl/c- \pcoBe-^\\v 3 \ > 

provided the no-ergosphere condition (]1.9p holds. The formula (I1.7P for the 
energy-momentum tensor turns into 

T a p= [ PaP ^(E,L)^=£l=, (2.2) 

where p a has to be expressed via (|2.ip . Here we first express the four 
dimensional integral in (jl.7p in terms of (v°, . . . ,v 3 ), replace the integration 
variable v° by m, and then use the fact that f is S distributed with respect 
to rn. In what follows we view / as a function on the mass shell. 

In order to turn on or off angular momentum we introduce a second 
parameter A £ 1, and in order to obtain the correct Newtonian limit for 
7 = we adjust the ansatz for / as follows: 

f = 4>{E-lh)i>{\L). (2.3) 
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The important point here is that 

e yu \/l + o'lfl 2 — 1 1 

E — I/7 = ¥ -Ivl 2 + 1/ + ujL as 7 -> 0; (2.4) 

7 2 

see (02) below. For 7 = this limit is to replace the argument of in (12. 3p . 
We specify the conditions on the functions and ?/;. 
Conditions on and 0>. 

(01) E C 1 (M) and there exists Eq > such that 0(??) = for 77 > £"o and 
0(77) > for 77 < Eq. 

(02) The ansatz f(x,v) = (^M 2 + U(x)), x,v £ M 3 , leads to a compactly 
supported steady state /at of the Vlasov-Poisson system, i.e., there 
exists a solution U = Ujy £ C 2 (]R 3 ) of the semilinear Poisson equation 

AU = 4tt Pn = 4vr J Q|u| 2 + J7J dv, 17(0) = 0, 

Un(x) = Un(\x\) is spherically symmetric, and the support of £ 
C 2 (1R 3 ) is the closed ball7J HiV (0) where U N {R N ) = E and U N (r) < E 
for < r < .Rat, C7v(r) > E for r > itV- 

(03) We have 

6 + 4irr 2 a N (r) > 0, rG [0,oo[, 

where 

a N (r) := / 0'f J |v| 2 + £fo(r)) du. 

(-0) -0 £ C C °°(M 2 ) is compactly supported, V > 0, <9l^(A,0) = for A G R, 
and -0(0, L) = 1 on an open neighborhood of the set 

{L = L(x,v) I (x,v) £ supp/at}. 

For the Newtonian steady state 

lim Un{x) = LTv(oo) > Eq. 

\x\— >oo 

The normalization condition C7v(0) = instead of Un(oo) = is uncon- 
ventional from the physics point of view, but it has technical advantages. 
Examples for ansatz functions which satisfy (01) and (02) are found in 
[6j[T3], the most well-known ones being the polytropes 

0(E) := (Eq - E) k + 
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for 1 < k < 7/2; here Eq > and (•)+ denotes the positive part. In m 
Sect. 7] it is shown that (4>3) holds for a subclass of the polytropes. 

We recall that the metric (jl.ip was written in terms of the axial coordi- 
nates p £ [0, oof, z G R, (p £ [0, 27r]. In what follows we shall also use the 
corresponding Cartesian coordinates 

x = (pcos f, psirnp, z) £ R 3 , 

and by abuse of notation we write u(p, z) = u(x) etc. It should be noted 
that tensor indices always refer to the spacetime coordinates t, p, z, <p. In 
Section [7] we collect the relevant information on the relation between regu- 
larity properties of axially symmetric functions expressed in these different 
variables. We can now give a detailed formulation of our main result. 

Theorem 2.1 There exist 5 > and a two parameter family 

(^7,A,-B 7 ,A,A i 7,A,^7,A)( 7 ,A)g[0,5[x]-(5,<5[ C C 2 (R 3 ) 4 
with the following properties: 

(i) (^0,05-80,05/^0,05^0,0) = (t^Vj 1)0,0) where Un is the potential of the 
Newtonian steady state specified in (4>2). 

(ii) If for 7 > a distribution function is defined by Eqn. i2.3\) and a 
Lorentz metric by M.l}) with c = 1/y/y then we obtain a solution of the 
Einstein- Vlasov system li 1.5\) . $1.6\) . {1. 7| ) which satisfies the boundary 
condition hi. 3^ and is asymptotically flat. For A / this solution is 
not spherically symmetric, and for appropriate choices of tp its total 
angular momentum is non-zero. 

(Hi) If for 7 = a distribution function is defined by Eqn. 12. 3\) . observing 
\2.1$ , this yields a steady state of the Vlasov-Poisson system with grav- 
itational potential z^o,a which is not spherically symmetric for A 7^ 0. 

(iv) In all cases the matter distribution is compactly supported both in phase 
space and in space. 

Remarks. 

(a) The smallness restriction to A implies that the solutions obtained are 
close to being spherically symmetric, and that their total angular mo- 
mentum is small. 
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(b) The smallness restriction to 7 = 1/c 2 is undesired because c is, in a 
given set of units, a definite number. However, if (f,u,B,pL,oj) is a 
solution for some choice of c G]0,oo[ then the rescaling 

Rp,z,p x ,p 2 ,p z ) = c~ i f(cp,cz,cp 1 ,cp 2 ,p 3 ), 
v{p,z) = c~ 2 v(cp,cz), 
B(p,z) = B(cp,cz), 
p(p,z) = p(cp,cz), 

u(p, z) = Uj(cp, CZ), 

yields a solution of the Einstein- Vlasov system with c = 1; notice that 
this rescaling turns (|1.9p into the corresponding condition with c = 1. 



(c) The metric functions v and p do not satisfy the boundary conditions 
P, but 

lim u(p,z)=u oc , lim p(p, z) = -v^/c 2 , (2.5) 

\(p,z)\-H)0 \(p,z)\^oo 

see Proposition 13.21 If we now abuse notation and redefine v = v — v^, 
p = p + Voo/c 2 and u = uje~ U(X / c then the original condition (|1.2p is 
restored and the metric (jl.ip takes the form 

ds 2 = - c ^/c\2 dt 2 

+ 4 ( e 2 »dp 2 + e 2 »dz 2 + p 2 B 2 e- 2 »/ c2 (d<p - LOddtf 



with constants c\ , C2 > which simply amounts to a choice of different 
units of time and space. By general covariance of the Einstein- Vlasov 
system (|1.5p . (jl.6p . (|1.7p the equations still hold. 

(d) In the course of the proof of the theorem additional regularity prop- 
erties and specific rates at which the boundary values at infinity are 
approached will emerge. 

We will transform the problem of finding the desired solutions into the 
problem of finding zeros of a suitably defined operator. The Newtonian 
steady state specified in (02) will be a zero of this operator for 7 = A = 0, 
and the implicit function theorem will yield our result. 

The Einstein field equations are over deter mined, and we need to identify 
a suitable subset of (a combination of) these equations which suffice to 
determine u,B,p,u, and which are such that at the end of the day all the 
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field equations hold once the reduced system is solved. To do so we introduce 
the auxiliary metric function 

£ = + fj,. 

Let A and V denote the Cartesian Laplace and gradient operator respec- 
tively. Taking suitable combinations of the field equations one finds that 



Au + — • Vu - -p 2 B 2 e^\Vu\ 2 = 4tt 
B T 11 



7 2 e ( 25 -4 7 „) (Too + 2a ;To 3 ) 



+7(Tn + T 22 ) + e 2 « f + J 2 u 2 e-^) T, 



33 



(2.6) 



AB + — • VB = 8irj 2 B (T xl + T 22 ) 
P 



Auj 



2 Yl + 3 ^ _ 4 7 V.) • V. = (T 03 + ^T 33 ) 



1 + 



a„5 



B 



5 P £ - p—^-ozi 

d p (p%B) - ^d zz B + 7 2 p {(d p uf - (d z vf) 



2pB 



1 + p- 



d n B 



B 



- 7 p 3 B 2 e- 4 ^ {(d p co) 2 - {d z u) 2 ) , 
a t , d * B a t 

dzi + p—{f d p{. 

dp{P l zB) + 2 1 2 pd p ud z v + ^BV 4 "^^. 
r> 2 



(2.7) 
(2.8) 



(2.9) 



(2.10) 



We write 



B = 1 + b. 



By taking a suitable combination of (|2.9p and (|2.10p we obtain equations 
which contain only d p £ or <9 2 £ respectively, and we chose the former. In 
the above equations the terms T a p are functions of the unknown quantities 
z/, b, u), £ = 7^+/i for which we therefore have obtained the following reduced 
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system of equations; throughout, = B, £, u), p; 7, A): 

Av = 4vr $00 + 7*11 + 27w$o3 + e 2? {l~^2 + 7 2 w 2 e~ 47 ^ $ 

_I V &.V^ + ^ 2 £V 4 ^|Vu,| 2 , 
2 

Ab + — • V6 = 8tt7 2 5$ii, 



33 



(2.11) 
(2.12) 



((i + b + pd p b) 2 + ( P d z b) 2 ) d P i 

= P d z b (d z b + P d zp b + irfpBdpvdzv + ^ 7 p 3 B 3 e~^d p ojd z u?j 

+(! + & + pd p b) ^(dppb + -b- d zz b) + 7 2 pB {{dpuf - (d z u) 2 ) 



1P 3 B 3 e~^ {(dpu) 2 - (d z cof 



(2.13) 



Aco + 2— • Vw 



16vr 



e 2? (7^03 + 7 2 w$33) - (^3^ - 4 7 Vz^ • Vw. 



(2.14) 



We supplement this system with the boundary condition (|1.3p . which in 
terms of the new unknowns reads 



£(0,z) = In (1 + 6(0, z)). 



(2.15) 



It remains to determine precisely the dependence of the functions Q a p on 
the unknown quantities u, b,£,u). The corresponding computation uses the 
new integration variables 



e 7 Vl + 7M 2 -l 



1 



We also introduce the notation 



m = 771(77, B, v, 7) = Be 



1 



, s = Be-^v 3 . 



/ e -2 7 f(i + 1V ) 2 - 1 

7 



l = l(s,B,v, 1 ) = -\ e^\ll + 1 



12 



and we obtain 

(DooO, B, e, u, p; 7, A) = 7 V«- 4 ^T 00 (2.16) 

27T Z" 00 f m 

-e 2 £~ A ' yu / / (1 + 777 + j psu}) 2 4>(rj + psu)ip(X, ps) ds drj 

J (e"i v -l)/7 J-m 



00 



= ^: e 2€-4 7 , / / (1 + 7??) 2 (/)(??)v , (A5 pa) d7? dSj 

$11 (i/,5,e,a;,p; 7, A) = Tn + T 22 (2.17) 

9 _ /'OO /*7Ti. 



„-e 4 / (m — s ) 4>(rj + psuj)i(j(\, ps)ds dr/ 

J{e'i v -l)h J-m 



B 3 



00 POO 



e 2 ^ I I (m 2 (i] — poos, B,u,j) — s 2 ) (j){rj)ip{X, ps)drj ds, 

00 J l+puis 



^ 33 (u,B,i,u,p; 1 ,X) = e 2 ^T 33 (2.18) 

2yj-p2 roc i-m 

= — — e 2 ^ / / s 2 <f)(7] + psuj)tjj(X, ps) ds dr\ 

B J( e yv_iyyJ_ m 

= B e2 l / s2 ^M\ ps) d V ds, (2.19) 

^ J— OO J l+pU)S 



<5> 03 {v, B, e, w, p; 7, A) = 7e 2? T 03 (2.20) 



2vrp ^ 



B 



/ s (1 + 777 + 7psw) 0(77 + psu))ip(X, ps) ds drj 



27TP p2( 

B 



(eT"-l)/7 J-m 
00 />oo 

e^* / / s (1 + 777) (j>(rj)tp(X, ps) drj ds; 

00 J l+pios 



we recall that T\\ = T 22 . In the course of the proof we will benefit from 
both of these two different representations of the matter terms. 

We now define the function spaces in which we will obtain the solutions 
of the system (|2.1ip - (|2.15p . By abuse of notation we write axially symmetric 
functions as functions of x G M 3 or alternatively of p > 0, z G R. We fix 
< a < 1/2, < (3 < 1 and consider the Banach spaces 

Xi := G C 2 ' Q (R 3 ) I u(x) = u(p,z) and \\u\\ Xl < 00}, 
X 2 := \b G C 3 ' a (R 3 ) I b(x) = b(p,z) and ||6||^ a < 00}, 
X 3 := {e G C 1 ' a (Z Ji ) I = 2) and < 00}, 

Af 4 := jw G C 2,a (]R 3 ) I uj(x) = uo(p, z) and IMI^ < 00 j, 
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where 

Z R := {x e R 3 | p < R) 

is the cylinder of radius R > 0, the quantity R being defined in (|2.21|) below, 
and the norms are defined by 

lkllo>.-(iii8) + ||(l + W) 1+ ^Vi/|| 00 , 

l|6|l03,a(RS) + ||(l + N) 3 V&|| 00 , 

II^Hc 1 ."^)' 

||w|| C 2,« (R 3) + ||(1 + kl) 4 Va;|| O0 , 
and 

X :=^x^ 2 x^x4 \{y,b&w)\ x := \\v\\ Xl + \\b\\ X2 + U\\ Xa + M\ Xi - 

Here || • || denotes the L°°-norm, functions in C k ' a (Q) have by definition 
continuous derivatives up to order k and all their highest order derivatives 
are Holder continuous with exponent a, and 

ll ll II rv ll \D°g(x)-D°g(y)\ 



\\ u \\x 2 
U\\x 3 



M\x 4 



|<x|</c \a\=k 



where D a denotes the derivative corresponding to a multi-index a £ Njj. It 
will be straightforward to extend £ to M 3 once a solution is obtained in the 
above space. 

The condition (4>2) on the Newtonian steady state implies that there 
exists R > Rjy > such that 

U N (r) > (E + U N (oo))/2 for all r > R. (2.21) 

If 

\\v - C/jvIL < \E Q - U N (oo)\/± and < 7 < 7o , 
with 70 > sufficiently small, depending on Eq and Un, then 

e ~/v(x) _ y 

> Eq for all Ixl > R. 

7 

Since L is bounded on the support of ip this implies that there exists some 
5q > such that for all {y, b, £, oj; 7, A) £ U the matter terms resulting from 
(|2.16p - (|2.18p are compactly supported in Br(0), where 

U :={(i/,6,e,w; 7 ,A) £ * x [0, <5 [x] - <5 , <5 [| 

- (U N ,0,0,0)\\ X < S }. (2.22) 
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In addition we require that 5q > is sufficiently small so that B = 1+6 > 1/2 
for all elements in U and the factor in front of d p £ in (|2,13p is larger than 
1/2. 

Remark. If we want to make sure that the no-ergosphere condition (jl.9|) 
holds for the solutions we construct we redefine 

MU == MIc-crs) + 11(1 + kl) 3 w|loo + 11(1 + M) 4 VHL. (2-23) 

If So > is sufficiently small then v is close to the given Newtonian potential 
Un, B is close to 1, and due to the redefined norm pa; is bounded so that 
(|1.9j) holds if 7 is sufficiently small. 

Now we substitute an element (z/, b, £, oj; 7, A) G W into the matter terms 
defined in (|2.16p - (l2.20p . With the right hand sides obtained in this way the 
equations ()2. 1 1 j) (|2. 14[) can then be solved, observing the boundary condi- 
tion (|2,15p . In order to do so we need to be a little careful with the definition 
of axial symmetry, because we will rewrite different equations as equations 
on W 1 with suitable, different dimensions n > 3. 

We call a function u : W 1 — > IR axially symmetric if it is invariant under 
all rotations about the x n -axis, or equivalently, if there exists a function 
u : [0, oo[xIR — y M such that 

u(x) = u(p, z), where p = \J x\ + . . . + x^ l _ l and z = x n for x G M n . 

Of course we will identify u and u. For what follows it is important that we 
can view an axially symmetric function u as a function defined on any W 1 
with n > 3. In particular, we remark that 

n— 2 Vp 
A n u = d P pU H d p u + 9 z2 n = An + (n — 3) — • Vu, 

where the left hand side refers to Cartesian coordinates on W 1 and the right 
hand side to M 3 . In view of this relation, equation (12.14p . when rewritten in 
terms of p and z, takes the form (jl.4p . The latter equation can be handled 
directly in these variables, cf. Lemma 17.21 but it is much more efficient to 
observe that this equation is nothing but a semilinear Poisson equation on 
IR 5 for an axially symmetric function. Similarly, (|2.12p is nothing but a 
semilinear Poisson equation on ]R 4 . 

Equation ()2.11|) and the properly rewritten equations ()2. 12[) and (|2.14p 
are solved (in terms of the right hand sides which of course contain the 
unknowns) by the corresponding Newton potentials, and (|2. 13[) can simply 
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be integrated in p. We define the corresponding solution operators by 



G 1 (u,b,^,uj;j,X)(x) :-- 



l 

+ — 

47T 



(y) 



dy 



G 2 (uA^lA)(x) :=-- / M 2 {y) T -^-^., 



, iEi 3 , (2.24) 
(2.25) 



|z-y| : 

rl 



G 3 (v,b,£,u>;j,\)(p,z) := hi (1 + 6(0,2)) + / g(s,z)ds, < p < R, 

Jo 



G A (v,b,£,u);j,X)(x) := —z [ 



Vb 

— ■ Vu + 47Vz/ • Vw - M 3 
B 



(2.26) 
dy 



x — y\ 3 ' 



(2.27) 



Here we put 



Mi : = ( $oo + 7$n + 270;$ 



n:; 



+ 7 



1 



p 2 B 2 

M 2 :=7 2 S$ 11 (^S,C,^, / 9; 7 ,A), 
16vr , , 2 



+ 7 Ve~ 47 " j $33 J (f , 5, £, w, p; 7, A) , (2.28) 

(2.29) 



3 



p 2 B 2 



(7^03 + rw$33)(v, B, C w, p; 7, A), 



(2.30) 



and 



g: = ((l + b + pdp b) 2 + (pd z bf) 1 

x P d z b(d z b + P d zp b + 2 1 2 P Bd p vd z v + i 7 p 3 J B 3 e - 47,/ d /3 wd z 

+ (1 + 6 + pd p 6) ( |( V + j^fe - a zz 6) + 7 2 pB ((c» 2 - {d z vf) 

- 1 p^B^e-^{{d p u J ) 2 - (d z uf) 



Finally we define the mapping to which we are going to apply the implicit 
function theorem as 

T :U -> X, (v, b,£,u; 7, A) h-> (1/, b,£,u) - (G\, G 2 , G 3 , Gi)(u, b, £, u>; 7, A). 
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In the next section we obtain the solutions in Theorem 12.11 as a two pa- 
rameter family of zeros of this mapping. It should be noted that functions 
resulting from the operators Gi will all be axially symmetric so that they 
can all be viewed as functions on M 3 even if they are at first defined on 
different M n 's. 

3 Outline of the proof 

We check in a number of steps that the mapping T satisfies the conditions 
for applying the implicit function theorem. Some of these steps turn out to 
be identical, or almost identical, to the corresponding steps in [4J in which 
cases the details will be left out. 

Step 1. 

We need to check that the mapping T is well defined. Since in this step the 
presence of to and also the somewhat different set-up for the space X causes 
some differences compared with the analysis in (3], we deal with this issue 
in some detail in Section [5l 

Step 2. 

The next step is to see that 

J^tv, 0,0,0; 0,0) = 0. 

This is due to the fact that for 7 = A = the choice b = £ = to = trivially 
satisfies (I2T2D . (I2T3D . (I2T4D . while (|2TTjl reduces to 

Au = 4tt<I>ooO, 1,0,0; 0,0) 

with 

$ooM, 0,0; 0,0) = 4tt^ HV)V2(V -v)drj = J <f> Qm* + J\ dv; 

notice that 6 = implies that B = 1. By (<j)2), v = Un is a solution of 
this equation, and the fact that Um £ X\ is part of what was shown in the 
previous step. Notice further that 7p(0, L) = 1 on the support of /at so that 
this factor, which vanishes for large L and formally makes the ansatz depend 
on L also in the Newtonian case, does not affect the Newtonian steady state 
at all. 

Step 3. 

The mapping T is continuous, and continuously Frechet differentiable with 
respect to (u, b, £, u). Since the new element to does not affect the proof for 
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the static case in any essential way we refer to [U Sect. 5] for the details of 
this step. 

Step 4. 

We have to show that the Frechet derivative 

C := DT{U N , 0, 0, 0; 0, 0) : X -»• X 
is one-to-one and onto. Indeed, 

C(5u, 5b, 5£, 5lo) = (5v - L x {8u) - L 2 (5b, 5£), 5b, 8£ - L 3 {5b), Suj 
where 

L\(8v)(x) ■= - ( 1 — - — 1 - |-r J a N (y)5u(y)dy, 
JM.3\\x-y\ \y\J 

L 2 (5b,50(x) :=±- [ V5b(y)-VU N (y) 



4vr J R 3 \x-y 
1 1 



\x-y\ \y\ 



PN{y)5£,(y) dy, 



L 3 (5b)(x) := 5b(0, z ) + \^ s (dppM + - p d P 8b ~ 9zzSb) (s, z) ds, 

< p < R, 

with oat as defined in (4>3). To see that C is one-to-one let L{5v, 5b, 5£, Su) = 
0. The second and the last component of this identity imply that 5b = = 
5uj, and hence also 5^ = by the third component. It therefore remains to 
show that 5v = is the only solution of the equation 5v = L\{5u), i.e., of 
the equation 

A5v = Aira N 5v, 5u(0) = 0, 

in the space X\. Under the assumption on ajy stated in (03) this was 
established in [U Sect. 7]. 

To see that C is onto let (g\, g 2 , 53,34) £ ^ be given. We need to 
verify that there exists [5v, 5b, 5£, Su) G X such that £(5u, 5b, S£, 5oj) = 
(<7i, g 2 , g 3 , 34). The second and fourth components of this equation simply 
say that 5b = g 2 and that 5u) = g&. Now 5b £ X 2 implies that L 3 {5b) E X% 
by Lemma [TTTT b) . Hence we can set <5£ = 53 + L%{5b) to satisfy the third 
component of the 'onto' equation, and it remains to show that 

Su - Li(5u) =gi+ L 2 (5b, 8£) 
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has a solution 8v € X\. Firstly, L2(6b, 5£) 6 X\. The assertion therefore fol- 
lows from the fact that L\ : X\ — > X\ is compact. We refer to [H Lemma 6.2] 
for the proof of this property. It is at this point that we use the fact that 
in X\ the decay assumption is weaker than what we actually get for Gi and 
that < a < 1/2: L\ gains some Holder regularity which is the source for 
the compactness. 

In view of the steps above we can now apply the implicit function theo- 
rem, cf. Thm. 15.1], to the mapping T :IA —> X . 

Theorem 3.1 There exists 61,62 £]0, 6[ and a unique, continuous solution 
map 

5: [0,$i[x] B S , 2 (U N , 0,0,0) CX 
such that 5(0,0) = (U N , 0,0,0) and 

J"(5(7, A); 7, A) =0 for all (7, A) e [0, <*i[x] - 6 U 6i[. 

The definition of J- implies that for any (7, A) the functions (y, b,£,oj) = 
5(7, A) solve the equations (|2.1ip - (|2.14p where the last equation deserves 
some explanation. By construction, 10 satisfies the equation Au = q on M 5 
where q is an abbreviation for the right hand side of (|2.14p . Both oj and 
q are axially symmetric, cf. Lemma 15.11 Hence uj and q can be viewed as 
functions of p and z, and as such they satisfy (jl.4p which in turn implies 
that as functions on IR 3 they satisfy (12. 14|) . An analogous argument applies 
to (gUP - 

If / is defined by ([23]) then the equations (f2TT3|) hold with 

the induced energy momentum tensor. We can extend £ to the whole space 
using the solution operator G3 for all x G M 3 . Also, the boundary condition 
(11.31) on the axis of symmetry is satisfied: 

C(0,z) = G 3 (v, b,0(0, z) =ln(l + b(0, z)) = In B(0, z); 

recall that £ = 7^ + \i. The solutions are asymptotically flat in view of 
Remark (c) given after the formulation of Theorem I2.lt a l so see Proposi- 
tion IO 

The solutions will in general have non-zero total angular momentum for 
7^0. To prove this assume that u> = 0. Then 

27T p f°° f' n 

T 3 = 5-/ (1 + 777)0(77) / stp(\,ps)dsdr). 

It is easy to see that there are functions ip satisfying the condition (-0) such 
that this integral is non-zero in contradiction to (I2.14p and u = 0, e.g., 
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choose tjj such that on the support of tp, tp(\,L) > ip(\, —L) for L > and 
A / 0. 

For 7 = we conclude first that b = 0, cf. (|2.12|) and the G2-part of the 
solution operator, respectively. Then the Cr3-part implies that £ = and 
the G^-part yields u) = 0, so that the solution reduces to (z/, 0,0,0), where 
v solves 

A*/ = 47r$ooM,0,0,p;0,A). 

Since 

roo ry/2{ri-v) 

*oo(f,l,0,0,p;0,A) =4tt / <p{r, + psu) if>(\, ps) ds drj 

Jv JO 

coincides with the spatial density induced by the ansatz f)2.3f) for the Newto- 
nian case, cf. [131 Lemma 2.1], part (iii) of Theorem 12. II is established. The 
resulting Newtonian steady state may or may not rotate, depending on the 
properties of ip\ see [131 Remark (b), p. 324]. 

The matter terms are compactly supported in view of of the discussion 
following ([221]). 

To complete the proof of Theorem 12.11 we must show that all the field 
equations are satisfied by the obtained metric (II. 1|) . The argument relies on 
the Bianchi identity V a G a ^ = which holds for the Einstein tensor induced 
by any (sufficiently regular) metric, and on the identity \/ a T a ^ = which is 
a direct consequence of the Vlasov equation (|1 .6j) ; V a denotes the covariant 
derivative corresponding to the metric (ll.ip . Due to the inclusion of the 
w-equation we cannot refer to the corresponding argument in [3], and we 
provide the details of the argument in Section [6l 

We conclude this outline of the proof of our main result by collecting 
some additional information on the solution which we obtain in the course 
of the proof and which shows that the solutions are asymptotically flat. 

Proposition 3.2 Let (u,b,^,cj) = 5(7, A) be any of the solutions obtained 
in Theorem\3J\and define p := i~vjc 2 and B := 1 + 6. Then £ G C 2 ' Q (IR 3 ), 
the limit := lim\ x \^. 00 u(x) exists, and for all a G Nq with \a\ < 2 and 
x G M 3 the following estimates hold: 

|L><>(z)-^oc)| <c(i + M)-( 1+H ), 

\D a (B- l)(x)| <C(l + |x|)- (2+H) , 
\D a £(x)\ <C Y (l + |x|)-( 2+H \ 
\D a u(x)\ < C(l + |x|r (3+H) . 
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In particular, the spacetime equipped with the metric is asymptotically 
flat in the sense that \2. 5\) and, after a trivial change of coordinates, also 
U.2\) holds; see Remark (c) after Theorem \2.1[ 

The proof follows easily from the decay estimates for various Newton poten- 
tials which we establish in Section [5l and from the fact that £ now satisfies 
the equations (|2.9p and (|2.1U|) ; see also [H Prop. 2.3]. In passing we remark 
that the asymptotic behavior stated in Proposition 13.21 agrees with what is 
given in [5]. The outline of the proof of our main result, Theorem 12.11 is 
now complete. 

4 Regularity of the matter terms 

In this section we investigate the regularity properties of the functions $oo> 
^33, and $03, and of the induced matter terms M\, M2, and M3 from 
(|238j) . (|239j) . and (l2~30l) . 

Lemma 4.1 Let <f> and tp satisfy the conditions (4>1), (4>2), and (ip). 

(a) The functions $00)^33; an d ^03 have derivatives with respect to 
£,v,uj,p, and B G]l/2,3/2[ up to order two, and these are continu- 
ous in v, £, B, lo, p, 7, A. For <J>ii the same is true with derivatives up 
to order three. 

(b) For (u, £, b, oj, 7, A) G U we have Mi, M 2 ,G C c 1,a (]R 3 ), M 3 G C C °' Q (IR 3 ), 
and Mi, M2, M3 are axially symmetric. 

Proof, (a) Differentiability with respect to £ is obvious to any order. Con- 
cerning differentiability with respect to v, uj, B, and p, the integrals in 
the formulae for expressed by means of / gain one derivative. Since 
4> G C 1 (IR) and tp G C~(M 2 ), the <3?ij have the desired regularity. For $n 
we have to observe the following fact. If we differentiate the integrand in 
the second form of (|2.17p with respect to one of the relevant parameters we 
obtain an expression which has the same structure as the in general. If 
we differentiate the integration boundary / + puis this gets substituted for n 
in the integrand and the term m 2 — s 2 vanishes. 

(b) By the choice of R in (I2.2ip and 5q in (I2.22D the matter terms which 
result by substituting an element from U into the $jj are compactly sup- 
ported. By the definitions of the spaces Xj the functions which are substi- 
tuted into &ij are axially symmetric and at least in C 1,a (]R 3 ). The expression 
$33 contains the factor p 2 so that the term (E^/p 2 , which is present in both 
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Mi and M 3 , lies in C7 1,a (M 3 ). Thus M 1 and M 2 belong to Cc' Q (M 3 ). In or- 
der to establish the assertion for M3 we only need to consider the expression 
N(p, z)j p where 

/OO /'OO 
/ s (1 + 777) <f>(r})ip(\, ps) dr\ ds. 
-OO J 1+pUJS 

We now think oiu,B,oj as functions of p G M and z € R which are even in p 
and lie in C 2 (R 2 ). Hence N G C 2 (]R 2 ), and iV is odd with respect topGl 
as can be seen by the change of variables sn-s. This easily implies that 
N/p is in C^{R 3 ), cf. 0J Lemma 3.2], and hence M 3 G C^M 3 ) C C c °' a (M 3 ) 
as claimed; notice also Lemma |7. II □ 
Remark. The additional regularity of <J>n is needed for the Frechet differ- 
entiability of G2] notice that G2 maps into C s,a . 



5 J 7 is well defined 

In this section we show that T(y, b, £, cj; 7, A) G Af for (v, b, £, cj; 7, A) G 
For the most part this is an assertion on certain Newton potentials. The 
following lemma collects the necessary information. 

Lemma 5.1 Let < a, S < 1, n G N mi/i n > 3, and 5 G C°' a (R n ) with 
\g(x)\ < C{l + \x\)~ n - 5 , x G R n . 

Define 

U(x):=- I , 9{V l 2 dy, xGR n 

T/ien £7 G C 2 ' Q (R n ), and for any a G Nft |cr| < 2, 

\D a U(x)\ < C(l + |x|) 2 - n -l CT l, x G R n . 
If g is axially symmetric, then so is U . 

Proof. Since g G L 1 n L°°(M n ) and Holder continuous, £7 G C 2 ' a (M. n ) with 

VC/(x) = (n-2) / 

J \x-y\ n 

and 

d Xl d Xj U(x) = (n- 2)a n 5ijg(x) 

+ (» - 2 ) / fty (^Rf) 0K») " <fo 
j|x-y|<d \F — 2/1 / 

+ (n " 2) / , ^(Sl)^^ 
■/|as— i/|>d \l x 2/1 / 
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for i,j £ {1, . . . , ra}, x 6 W 1 , and d > 0, cf. [8]; a n denotes the volume of 
the unit ball in IR n . We consider the decay of the gradient first: 



|W(x)| < (n-2) 
< C 



\g{y)\ 



x-y\<\x\/2 \ x ~ V\ 
\x-y\<\x\/2 



L_ d y + („ _ 2) 

dy 



( 1 + \y\r n ' 5 w ^z- 1 +c\\ 9 \\i\x 



x-y\>\x\/2 
1-n 



l-n 



<C(l + N/2)— s [ T^^I + C \ X \ 

J\x-y\<\x\/2 \ x ~ U\ 

< C(l + \x\)- n ~ 5 \x\/2 + C\x\ l ~ n < Clxl 1 "™. 



The decay for U follows in the same way. For the second order derivatives 
we observe that the first term in the formula above decays as required by 
assumption on g. We choose d = (1 + |x|)~ n / Q . Then 



f d ( xi-yi \ 

J\x-y\<d Xj y\ n ) 



ia{y) - g{x)) dy 



< C||#|| c o, Q ( ffi 3) 

= c(\ + \x\r n 



dy 



\x-y\<d 



\x - y] 1 



C 



~ a - l dr = Cd a 



In order to estimate the remaining term we consider x £ W 1 with |x| suffi- 
ciently large: \x\/2 > (1 + \x\)~ n / a . Then 



I 







X% Vi 



y \>d 3 \\x - y\ n 



< C 

< c 



9{y) dy 

\g{y)\ dy 

d<\x- y \<\x\/2 \x - y\ n 

■ , , + \x vV 

d<\x-y\<\x\/2 \X — y\ 



+ C 



-n—S 



x-y\>\x\/2 

dy 



+ C\\g\\i\x\ 



< C(l + \x\/2)- n ~ 5 [ 

Jd<\x 



dy 



d<\x— y\< 



l*|/2 \x-y\ 7 



+ C\x\~ n 



< C(l + \x\/2)- n ~ s In + \x\) n/a J + C\x\~ n < C\x\- n . 

Assume now that g is axially symmetric, i.e., the function is invariant under 
rotations about the Then the equation AU = n(n — 2)a n g, which is 

satisfied by U on M n , is invariant under these rotations. Since U vanishes at 
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infinity, it is the unique solution of this equation and hence axially symmetric 
as well. □ 



Lemma 5.2 Let {y, b, £, u; 7, A) G U and let G\ = G\{y, b, £, u); 7, A) 6e as 
de/ined in fOfl }. Taen G x G X x 

Proof. The source term M\ of the first part of G\ is in Cc' a (]R 3 ) by 
Lemma 14.11 (b). The definitions of Xi, X2, and X4 imply that the source 
term 

Iv&.Vi/- V^ 2 e- 4 ^|Vu;| 2 

of the second part of G± lies in C°' Q (]R 3 ) and decays like (1 + |x|) _3_/3 . 
Moreover, both source terms are axially symmetric. Hence Lemma 15 . 1 1 with 
n = 3 implies the assertion. □ 
We notice that in the proof above we did not need to use the full available 
regularity. 

Lemma 5.3 Let (v, b, £, u; 7, A) € U and let G4 = G^v, b, £, w; 7, A) 6e as 
de/med in fP7| J. T/ien G4 £ A4. 

Proof. We can view the source term 

V6 

q := M 3 - 3— • Vw + 47VZ; • Vw 
B 

as an axially symmetric function both on IR 3 and on M 5 . By Lemma [4. II (b) 
and the definitions of X\, X2, and X4, q G C 0,a (IR 3 ) and hence also q G 
C 0,a (lR 5 ). Moreover, the compact support of M3 and the decay estimates 
for Vz>, V6, and Vw in the corresponding spaces imply that q decays like 

(i + \x\)- 5 -? ~ (i+ p +|*|)-5-0 

where x G M 3 or in M 5 has axial coordinates p and z. Lemma ED with n = 5 
implies the assertion. The latter is true with the original definition of the 
norm in X4 as well as with (I2.23p . □ 

Lemma 5.4 Let (u, b, £, oj; 7, A) G U and let G2 = G2(v, 6, £, w ! 7) A) 6e as 
defined in 112. 25\) . Then G2 G <%2- 

Proof. The source term M2 lies in Cc' a (M 3 ), and since it is axially symmet- 
ric we can equally well view it as a function in Cc' a (M 4 ). Thus Lemma 15. II 
with n = 4 implies the assertion; notice that here we can throw one deriva- 
tive onto the source term so that G2 ends up in C 3,a . □ 

It remains to see that also G3 G X%, but the corresponding proof is 
identical to the one in [U Lemma 4.2 (c)]. 
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6 All field equations hold 



Lemma 6.1 Let (v,b,fj, = £ — 7^, oj) be one of the solutions obtained in 
Theorem \3. 11 Then the metric U.l\) together with f defined by \2. 3\) solve 
the full Einstein- Vlasov system U.5\) . hi. 6^ , and fl 1. r /\ ). 

Proof. For a metric of the form the components 00, 11, 22, 33, 03, 

and 12 of the field equations are nontrivial. We have so far obtained a 
solution v,B,£,uj of the reduced system ([2TL11 . (I2TT41 . (I2TT21) . (I2TT31) . We 
define E a p := G a p — 87rc _4 T aj g so that the Einstein field equations become 
E a p = 0. By the reduced system, 

E 00 + c 2 e- 2 ^ c2 \E n + E 22 ) + 2u J E 03 +^-^e^ c2 + c/) £33 = 0, (6.1) 

E 03 + uE 33 = 0, (6.2) 
E u + E 22 = 0, (6.3) 

1 + p^f\ (E n ~ E 22 ) + p^-E 12 = 0. (6.4) 

The Vlasov equation implies that V Q T Q ^ = 0, and V a G a ^ = due to 
the contracted Bianchi identity where V Q denotes the covariant derivative 
corresponding to the metric (jl.ip . Hence V a E a ^ = 0. We want to use 
these relations for /? = 1 and (3 = 2 together with (|6.ip - (|6.4p to show that 
E a /3 = 0. To do so we first rewrite the equations (|6.ip - (16,4p in terms of E a @ . 
Then (E21) and (EH turn into 



uE 00 _ £03 = Q) (6.5) 
E u + E 22 = Q (6.6) 
Using these to eliminate E 22 and E" 03 the equations (|6,ip and (16.4p become 
c 2 e Mc 2 _ ^bV) £ 00 + p 2 ^ 33 = 0, (6.7) 



L + P^W+/^ 12 = 0. (6. 
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The two Bianchi equations mentioned above can be written in the form 

d p ( P E^) + d z {pE™) + (Ad p p + {pE ^) + [Ad z p + *jf) (pE a ) 

+ p (r^o + 2^ 3 ) E 00 + pTl 3 E 33 = 0, (6.9) 

d p ( P E^) - d,{pE*) + + {pE") - (Ad z p + ^) ( P E^) 

+ p {T 2 + 2 W r§ 3 ) + ^2^33 = Qi (6 1Q) 

where (|6.5p and (|6.6p were used to eliminate -B 03 and E 22 . 

At this point there is a small subtlety concerning regularity. By definition 
of X% we have £ G C 1 '", and since this function satisfies (I2.13P also d p ^ £ 
C 1,a . Hence d p (pE u ) and d p (pE 12 ) exist classically, but a priori this need 
not be true for d z (pE u ) and d z (pE 12 ). However, approximating £ by smooth 
functions the corresponding Bianchi identities again hold, and passing to the 
limit, (|6.9p and (16.10P hold in the sense of distributions. With the possible 
exception of d z (pE 11 ) and d z (pE 12 ) all the terms in (|6.9p and (|6.10p are 
continuous so that these identities show that the latter terms are classical, 
continuous derivatives as well. 

We use (|6.7p and (16. 8p to eliminate E 00 and E n from the two Bianchi 
identities (|6.9p and (|6.10p . The resulting two equations contain only E 12 
and its first order derivatives and E 33 . Eliminating the latter finally yields 
the following first order partial differential equation for pE 12 : 

((->W+K^)>* 12 > 

- \ (l + P & -f) (p^) d z ( P E' 2 ) + c(p, z) (pE^ 2 ) = 0. 

Here c = c(p,z) is a continuous function on [0, 00 [xR the form of which is 
of no further interest, and the equation holds for p > 0. 

We recall that by definition (|2.22p of the set U the quantity 1 + pd p B / B 
is bounded away from zero so that any characteristic curve of the above 
equation must intersect the axis p = where pE 12 = e~ 2 ^(d z B / B — "yd z v — 
d z p) vanishes due to ([L3]) . Hence E 12 = on [0,oo[xR. By (jSH]) also 
E u = 0, and by flEU) the same is true for E . If we eliminate E 00 from 
(|6.9p using (I6.7P we find an equation of the form (1+ pd p B / ' B)E 33 = so that 
E 33 = 0. By dfTTJ) this implies that £ 00 = provided we choose the smallness 
parameter in the definition (|2.22|) of the set U such that the coefficient 
of E 00 in (]6.7p does not vanish which is possible due to the estimates in 
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Proposition 13.21 By (|6.5p we finally find that E 03 = 0, and hence all the 
Einstein equations E a p = hold. □ 

7 Appendix: On the regularity of axially symmet- 
ric functions, and a remark on the u equation 

We first collect a some remarks concerning the regularity of axially symmet- 
ric functions on W 1 , n > 3. 

Lemma 7.1 Let u : M. n — > R be axially symmetric and u{x) = u(p, z) where 
u : [0, oo[xR -> R. Let k G {1, 2, 3} and a g]0, 1[. 

faj u G C fc (R n ) <J=^> u G C fc ([0, oo[xR) and a// derivatives of u of order 
up to k which are of odd order in p vanish for p = 0. 

(b) u G C°' a (M n ) ^ C°' a ([0,oo[xR). 

The case n = 3 is already stated in [H Lemma 3.1], and the proof does not 
depend on the space dimension. 

We conclude this paper by pointing out that the cu equation (|1.4p can 
also be solved directly in the original variables p and z. However, as the 
authors had to learn, basing the analysis on the following result makes it 
much harder. 



Lemma 7.2 A solution to the equation ^l-4\ ) is given by 
u(p,z)= / / JC(p,z,p,z)q(p,z)dzdp, 



where 

^^--^{p) 2p~p J' 

and Q1/2 is a half-integer Legendre function of the second kind. 

Sketch of proof. As explained above, we can interpret q as an axisymmet- 
ric function on R 5 and rewrite (|1.4p as the Poisson equation on R 5 . Hence 
the solution can be represented as 

1 f l(y) j.. _ in.5 



u(x) = -—s / | prdy, x G 

8vH J r5 \x - y\ A 
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If we let x = (0, 0, 0, p, z) and introduce polar coordinates for the integration, 
two integrals can be carried out explicitly, and the fact that 
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